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ABSTRACT 

Theoretically expected natures of a supernova driven by a wind/jet are discussed. 
Approximate analytical formulations are derived to clarify basic physical processes 
involved in the wind/jet-driven explosions, and it is shown that the explosion prop- 
erties are characterized by the energy injection rate (Ei so ) and the mass injection 
rate (M- lso ). To explain observations of SN 1998bw associated with Gamma-Ray Burst 
(GRB) 980425, the following conditions are required: E iso M iso > 10 51 erg M© s~ 2 and 
-Eiso > 2 x 10 52 erg s _1 (if the wind Lorentz factor T w ~ 1) or E iso > 7 x 10 52 erg 
s" 1 (if T w > 1). In SN 1998bw, 56 Ni (~ O.4M ) is probably produced in the shocked 
stellar mantle, not in the wind. The expected natures of SNe, e.g., ejected 56 Ni masses 
and ejecta masses, vary depending on E[ so and M- mo . The sequence of the SN prop- 
erties from high E lso and Mj so to low E mo and M- m o is the following: SN 1998bw-like 
- intermediate case - low mass ejecta (< 1M©) where 56 Ni is from the wind - whole 
collapse. This diversity may explain the diversity of supernovae associated with GRBs. 
Our result can be used to constrain natures of the wind/jet, which is linked to the 
central engine of GRBs, by studying properties of the associated supernovae. 

Key words: gamma-ray: bursts - supernovae: general - supernovae: individual (SN 
1998bw) - nuclear reactions, nucleosynthesis, abundances. 



1 INTRODUCTION 

Gamma-Ray Bursts (GRBs) are energetic cosmological 
events, emitting > 10 51 ergs in 7-ray. A leading model for 
the central engine of GRBs is the formation of a black hole 
(BH) and an accretion disk, following the gravitational col- 
lapse of a massive star whose main sequence mass (M ms ) 
is at least as large as 25Mq (for reviews, see Woosley & 
Bloom 2006; Nomoto et al. 2007). A relativistic flow gen- 
erated by neutrino annihilation (Woosley 1993; MacFadyen 
& Woosley 1999) or magnetic activity (Brown et al. 2000; 
Proga et al. 2003) is proposed to trigger a GRB. 

A link between (a class of) GRBs and Type Ic super- 
novae (SNe Ic) has been established observationally. The 
most convincing cases for the supernovae associated with 
GRBs (hereafter GRB-SNe) have been provided by spectro- 
scopic detection of supernova features in an optical afterglow 
of a GRB or at the position consistent with a GRB. Three 



nearby GRB-SNe detected in this way are found to be sim- 
ilar to one another. The category includes GRB 980425/SN 
1998bw (the proto-typical GRB-SN; Galama et al. 1998), 
GRB 030329/SN 2003dh (Hjorth et al. 2003; Kawabata et 
al. 2003; Matheson et al. 2003; Stanek et al. 2003), and GRB 
031203/SN 20031w (Malesani et al. 2004; Thomsen et al. 
2004). Optical observations of these GRB-SNe are well ex- 
plained by an explosion of a carbon-oxygen (CO) star, which 
has evolved from a massive star (M ma ~ 4OM0) and has 
lost its H- and He-envelopes during the hydrostatic evolu- 
tionary phase (Iwamoto et al. 1998; Woosley et al. 1999; 
Nakamura et al. 2001a; Mazzali et al. 2003, 2006). The ki- 
netic energy (Ek) of the expansion is large, E51 = Ek/W 51 
ergs > 10 (note that E51 ~ 1 for canonical supernovae). 
They eject ~ 0.3 - O.7M of 56 Ni (which powers the SN 
luminosity by the decay chain 56 Ni — > Co — > Fe). Hereafter, 
the mass of 56 Ni is denoted by M( 56 Ni). Recently, another 
example of the association has been reported (Delia Valle 
et al. 2008; Soderberg et al. 2008), i.e., GRB 081007/SN Ic 
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2008hw, while the observed properties of this SN have not 
been modeled yet. 

Despite the similarity within the well studied cases men- 
tioned above, GRB-SNe do seem to have diverse properties. 
Peak magnitudes of so-called supernova bumps seen in GRB 
optical afterglows show diversity (Zeh, Klose, & Hartmann 
2004; Woosley & Bloom 2006), highlighted by sub-luminous 
(possible) SNe in GRBs 040924 and 041006 (Soderberg et 
al. 2006). A few GRBs show no evidence for the supernova 
bump (Hjorth et al. 2000; Price et al. 2003). Non-detection 
of SN features in two nearby GRBs 060505 and 060614 has 
been reported, placing the upper limit to brightness of pos- 
sible underlying SNe ~ 100 times fainter than SN 1998bw 
(Delia Valle et al. 2006; Fynbo et al. 2006; Gal- Yam et al. 
2006). 

In spite of the observational constraints, the explosion 
mechanism of GRBs and GRB-SNe are still unknown. In 
particular, how properties of the central engine are related 
to the bulk expansion of the stellar materials, observed as 
a supernova mainly in visual light, is still under debate. A 
possibility is that a supernova is induced by a disk wind 
generated by viscous heating (MacFadyen & Woosley 1999; 
Narayan, Piran, & Kumar 2001). 

Although there were numerical calculations for the 
wind/jet-driven explosions (Khokhlov et al. 1999; Mac- 
Fadyen, Woosley, & Heger 2001; Maeda & Nomoto 2003; 
Nagataki et al. 2003; Maeda 2004, Nagataki, Mizuta, & Sato 
2006; Tominaga 2007a; Tominaga et al. 2007b; and Tomi- 
naga 2009), it has not been clarified what fundamentally de- 
termines the properties of resulting SNe and in what ways. 
Also, the numerical investigations have been restricted in 
the parameter space. 

Aiming to overcome these problems, this study is com- 
plementary to the past numerical studies. Our goal in this 
paper is to express theoretically expected features of SNe 
resulting from the wind/jet-driven explosion, as a function 
of rates of the mass and energy (M w and E w , where the 
subscript "w" denotes "wind", or Mi so and -Ei so , referring to 
the isotropic equivalent values) generated and injected from 
the central system (i.e., a black hole plus a disk) into the 
surrounding stellar mantle (in this paper, stellar "mantle" 
refers to the stellar materials above the central remnant, 
i.e., the outermost layer of the Fe core, and the Si- and CO- 
layers) . 

Our strategy is the following. (1) We first clarify what 
are main ingredients of the wind/jet-driven SN explosion 
(§2). We develop a simplified description for the shock prop- 
agation and nucleosynthesis in the explosion (§3). We es- 
pecially focus on the production of 56 Ni, addressing how 
two proposed sites for the 56 Ni production, a shocked stellar 
mantle and a disk wind (MacFadyen & Woosley 1999; see 
Maeda & Nomoto 2003 for a review), can be distinguished. 
(2) We then compare the results with observations of SNe as- 
sociated with a GRB, in order to constrain Mi so and _Ei so in 
these SNe (§4). The required values for Mi so and -Bi so then 
should be regarded as conditions that any models for the 
central engine should satisfy. In other word, by constrain- 
ing Mj ao and i5i so , we aim to provide useful constraints in 
studying the properties of the central engine of GRBs. 



2 MODELS 

We consider a situation that a supernova explosion is driven 
by outflow of materials (e.g., a disk wind; MacFadyen & 
Woosley 1999) from the vicinity of a central remnant (likely 
a black hole). This energy input can be different from the 
relativistic jet producing a GRB; The wind can either be 
relativistic or non-relativistic at its injection from the central 
remnant into the surrounding stellar mantle. 

A schematic picture of the problem considered in this 
paper (and most of the past numerical studies) is shown 
in Figure 1. Throughout this paper, we adopt 25Mq and 
40M© progenitor models from Nomoto & Hashimoto (1988). 
Following the treatment of the past numerical studies, the 
wind/jet (we hereafter frequently call it simply a wind) is 
injected by hand at a certain radius (i? w in this paper). 

2.1 Properties of the wind/jet 

In a one- dimensional hydrodynamic problem, three inde- 
pendent variables must be specified (two for thermody- 
namic variables and one for a hydrodynamic variable) as 
initial/boundary conditions. For simplicity, we assume that 
the wind is dominated by the kinetic energy. This does not 
drastically defeat our results, since it is expected that the 
thermal energy deposited at the root of the wind, near the 
central remnant, is quickly converted to the kinetic energy 
well below i? w . Thanks to this simplification, only two inde- 
pendent variables at i? w are required to determine the hy- 
drodynamic evolution of the system. The choice of the inde- 
pendent variables can be arbitrary, thus we take the energy 
injection rate (i5 w ) and the mass injection rate (M w ). The 
Lorentz factor of the wind/jet at the injection (at i? w ) is ex- 
pressed as a function of these two independent variables, i.e., 
r w ~ l+i? w /(M w c 2 ), where c is the speed of light. Hereafter, 
quantities expressing the properties of the wind at injection 
(at i? w ) are denoted by the subscript " w" . Trivially, the wind 
is initially highly relativistic only if -E w ,5i/M Wj ^> 1000, 
where E m 51 = E w /(10 51 erg s^ 1 ) and M w Q = M W /(M© 

Our analysis is based on a one-dimensional radial flow; 
The collimation of the wind is taken into account with a 
geometrical factor /n(^ 1). fn relates the wind intrinsic 
properties and isotropic equivalents as 

Bw = fnEi BO , (1) 
Mw = /oMiso , (2) 
E v = fnEi BO , (3) 

where _E W is the intrinsic total energy (i.e., the injected en- 
ergy integrated over time), and the subscript "iso" refers to 
isotropic equivalent values, fn is a measure of the collima- 
tion angle of the wind at 7? w : When the wind is a spherically 
symmetric flow, fn = 1. The wind is more narrowly colli- 
mated at the injection for smaller fn- 

For the temporal history of the wind properties, we fo- 
cus on the situation that the wind is injected at i? w , with _B W 
and M w constant in time for t < t w , and then terminated 
at t = i w , so that 

E w = E„t„ (or equivalently, E iso = i5i so iw) ■ (4) 

Some of the following results, however, do not rely on this 
assumption. 
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Figure 1. A schematic picture of a wind/jet-driven SN explosion, (a) A situation before the launch of the outward shock wave. The 
wind/jet is assumed to be continuously injected from the central remnant with properties denoted by the subscript "w" at radius R w . 
When materials initially at M r in the presupernova progenitor fall to -R w , their properties are denoted by the subscript "m". At this 
point, the mass of the central remnant (the mass within i? w ) is M r . (b) A situation after the launch of the outward shock wave. The 
shock wave arrives at radius -R sw , and temperature behind the shock is denoted by T. 56 Ni is produced both in the shocked stellar 
materials between R sw and i? w and within the wind/jet itself. 



2.2 Processes involved 

2.2.1 Collapse to Explosion 

The first function we have to consider is the dynamical effect 
(§3.1) of the wind to the collapsing stellar mantle; the over- 
lying stellar mantle continues to collapse onto the central 
remnant, and the wind does not always have a sufficiently 
large momentum to overcome the ram pressure of the in- 
falling materials (Fig. la). The importance of the dynam- 
ical effect to determine the outcome of a wind/jet-driven 
supernova explosion was pointed out by Maeda (2004) and 
numerically examined by Maeda (2004) and Tominaga et al. 
(2007b). Some numerical calculations have included the dy- 
namical effect self-consistently (e.g., Maeda & Nomoto 2003; 
Maeda 2004; Tominaga et al. 2007b; Tominaga 2009), while 
others have not (e.g., Maeda et al. 2002; Nagataki et al. 
2003). 



difficult to follow the dynamics of the shock wave propagat- 
ing into the stellar mantle; Thus, it has not been yet clear 
how the amount of 56 Ni synthesized in the wind is connected 
to properties of the progenitor star and of the resulting su- 
pernova. 



2.2.3 Geometry of the ejecta 

Finally, the properties of the wind affect the shape of the 
supernova ejecta, as the shock wave propagates through the 
stellar mantle. This has been directly examined in the past 
numerical studies, but restricted in the parameter space. In 
this paper, we derive a simple estimate of the shape, in terms 
of the properties of the wind (§3.3). 



2.2.2 Production of ° 6 Ni 

Once the shock wave is launched, it propagates outward 
into the stellar mantle (Fig. lb). The kinetic energy is now 
converted to the thermal energy following the shock wave 
propagation. The temperature behind the shock wave is ini- 
tially so high that the nuclear burning converts the initial 
stellar composition (mostly oxygen and silicon) mainly to 
56 Ni (§3.2.1). The temperature decreases as the shock wave 
moves outward; once the temperatures drops below ~ 5 x 10 9 
K, then the efficiency of the production of 56 Ni decreases 
rapidly. 

At the same time, the temperature of the injected wind 
itself may also be sufficiently high so that a fraction of this 
material may be converted to 56 Ni (§3.2.2). In most of the 
past numerical studies, the production of 56 Ni in shocked 
stellar mantle has been investigated in detail. The ejection 
of 56 Ni in the wind has been examined by a different kind 
of numerical simulations involving the innermost part of the 
collapsing star (e.g., MacFadyen & Woosley 1999; Nagataki 
et al. 2007). However, in such simulations, it is practically 



3 WIND/JET-DRIVEN EXPLOSIONS 
3.1 Collapse to Explosion 

In this section, we examine dynamical effect of the wind 
on the collapsing stellar mantle (Fig. la). We follow analy- 
sis similar to that given by Fryer & Meszaros (2003) for a 
standard neutrino-driven delayed explosion for (canonical) 
supernovae. 

We evaluate outcome of the interaction between the 
wind and the infalling materials at radius i? w ~ 8 x 10 7 
cm (for M ms = 40M Q ) and 1.2 x 10 8 cm (M ms = 25M ), 
where the presupernova enclosed mass (M r ) is 1.4M . As 
time goes by, materials initially at larger M r collapses to 
the radius R w , and add to the mass of the central rem- 
nant (Mbh). During a whole period before the launch of the 
shock, the temporal evolution of the system is thus speci- 
fied by the central remnant mass Mbh which monotonically 
increases as a function of time (i.e., larger Mbh for later 
time). 

If the trajectory of the infalling materials is that of free 
fall, the density (p m ) and velocity (w m ) of the infalling ma- 
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Figure 2. Density Structure (left vertical axis label) following 
the gravitational collapse, at Mbh = 2.4Mg (thick-solid line), 
3AMq (medium-solid), and 13Mq (thin-solid). Also shown is the 
postshock temperature (Tg = T/10 9 K; right vertical axis label) 
as a function of the shock radius, for Mbh = 2AMq, Ei SOi 5i = 30, 
and E lso = 10 (dashed line) or 1 (dotted line). The progenitor 
model (M ms = 40M Q ) is from Nomoto & Hashimoto (1988). 
The position within which the postshock temperature is above 
5 X 10 9 K (horizontal thin-dashed line), for M B h = 2.4M , is 
marked by the open circle (Bi so ,5i = 10) and by the open square 
(-Eiso,5l = 1). 



terials at R v are written as 



Pprcsn 



( R prcsn \ 

V R„ J 



3/2 



2GM r 



( Rv 



1/2 



(5) 



(6) 



The subscript "presn" is used for the pre-collapse initial 
values for the material at the mass coordinate M r . The 
subscript "m" refers to quantities of the same material (at 
M r ) when it collapses to i? w (Figure la). Figure 2 shows 
snap shots of the density structure for different Mbh for 
M ms = 40M Q . 

The outward shock wave is launched if the ram pressure 
of the wind (P w ) overcomes that of the infalling materials 
(P m ), i.e., 



fw pw^w ^ Pm 



pmV„ 



Pm9- 



GMbh 



(7) 



Here M B h is the mass of the BH when the material at M r 
(— M B h) collapses to i? w . Although a numerical constant 
g = 2 in the above estimate, we find g = 1/2 yields a bet- 
ter representation of a set of numerical simulations (Maeda 
2004; Tominaga et al. 2007b). We assume g = 1/2 through- 
out the paper. 

Equation (7) can be expressed in terms of E„ and M w , 
which are related to other properties of the wind/jet as fol- 
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Figure 3. The requirement for Ei BO to initiate the explosion as a 
function of M B h- For T w ~ 1, the case with M iso = O.IMq s _1 
is shown for presentation (note that the product .Ej so Mj so should 
exceed the value corresponding P m if T w ~ 1; equation 10). Three 
horizontal dashed lines at -Ei SOj 5i = 50, 10, and 7 are shown 
for illustrating purpose (arrows indicate the position where the 
outward shock wave is launched); Let us assume M ms = 40Mq, 
r w ~ 1 and M iso = O.IMq s _1 . If E iBOt5 i = 50 or 10, the outward 
shock wave is launched at the intersections between the Si- and 
CO-layers at M B h ~ 2.4M . On the other hand, if Ei so ,5i = 7, 
the momentum can not exceed the infalling ram pressure until 
nearly the whole CO star collapses to the central remnant at 
M BH ~ 13M Q . 



lows. 

M w 



(r w — l)p w f w c , and 

^4 w r w p w l> w , 



(8) 
(9) 



where A w is the area subtended by the wind. The asymp- 
totic expressions for the momentum balance are derived by 
substituting equations (8) and (9) into equation (7); 



S w M n 



fnV -EisoMiso > -^=A w p ln g G ^ BH 
V2 K w 



for r w 



for T w > 1 . 



GM B 



(10) 



(11) 



The RHS's of equations (10) and (11), except for A„ de- 
pending on fn , are completely determined by the progenitor 
structure, as a function of Mbh (Fig. 2). The requirement for 
E„ and M w in terms of _Ei so and Mi so is shown in Figure 3 
as a function of Mbh (i.e., as a function of time). Hereafter, 
we use a notation X n = X/10 n in CGS unit. 

If T w ~ 1, it is necessary that the product Ei so Mi so 
should be larger than a specific value given by equation (10), 
which corresponds to P m as a function of Mbh, in order to 
initiate the explosion. For T w 3> 1, -Eiso should exceed the 
value given by equation (11), and the requirement becomes 
independent from Mi so . In other ward, once the temporal 
evolution of Ei BO and Mi so is given, the outward shock wave 
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is launched when Ei so Mi so (or Ei so if T w 3> 1) is higher than 
the specific value representing P m (Fig. 3). 

Although P m overall decreases as a function of Mbh 
following the density decrease, jumps are seen at the edges 
of the characteristic hydrostatic burning layers. Assuming 
freefall, P m oc p P rcsn-R P rcsnMBH. Within each layer, p proS n 
drops slowly as a function of Mbh and J? pre sn, making Pm 
nearly constant as a function of Mbh- On the other hand, 
at the edg 6S, Ppresn decreases suddenly, leading to rapidly 
decreasing P m . As a result, local minimum values appear at 
the intersection of different layers (Fig. 3). It infers that the 
explosion is likely initiated at one of these intersections, but 
not within each layer. 

This behavior of P m provides an interesting implica- 
tion. For illustration purpose, let us take the case with 
M ms = 40Mq and r w ~ 1, assuming that Ei SO ,5i and 
Mi so are constant in time. If Eiso.siMso.o > 1, then the 
explosion is initiated below the Fe/Si interface (Mbh < 
2.4Mq). If 0.8 < -Biso,5iMi SOi < 1, the explosion position 
jumps to the Si/CO interface (Mbh ~ 3.4M ). Then, if 
Siso,5iMi so ,o < 0.8, almost whole CO layer collapses onto 
the central remnant. Thus, we conclude that a small dif- 
ference of the wind properties can lead to totally different 
outcome, with the critical value Eiso,5iMis O ,0 ~ 0.8 — 1 (for 
T w ~ 1) or E iso ,5i ~ 70 (for T w > 1). The critical value de- 
pends on M ms , and is smaller for a less massive progenitor 
(Fig. 3). 



3.2 Production of 56 Ni 

Once the ram pressure of the wind (P w ) overcomes the ram 
pressure of the infalling material (Pm), an outward shock 
wave is launched and explosive nucleosynthesis takes place. 
In this section, we examine production of 56 Ni as the shock 
wave passes through the stellar mantle (Fig. lb). We first 
discuss production of 56 Ni in a stellar mantle heated by the 
shock wave (§3.2.1), then we comment on production of 56 Ni 
within the wind/jet (§3.2.2). 



3.2.1 Shocked stellar mantle 

The shocked mantle can become the predominant site for the 
56 Ni synthesis, if the shock wave sweeps up a large amount 
of the stellar mantle, at least ~ O.IM©. This inevitably re- 
sults in a non-relativistic shock wave, even if the wind/jet 
at its emergence from the central system at R v is highly 
relativistic. Let us denote the average velocity and the mass 
of the expanding material swept up by the shock wave, by 
r shock and M shock , respectively, then 



100 



r sho ck ~ 1 + 6 x 10 



— 4 Ev, 



it 



M, 



shock,© 



(12) 



at time t (in second). In order to accelerate ~ O.IMq of 
the stellar mantle materials to r s hock = 100, the explosion 
energy is required to be higher than 10 55 ergs 

We can thus use non-relativistic approximation for the 
postshock temperature (T) as a function of the radius of 
the shock wave (i? S w). For a radiation-dominated fireball 
expanding into a uniform medium with the density p (Maeda 
& Nomoto 2003), 
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Figure 4. The isotropic mass of 56 Ni synthesized in the shocked 
stellar mantle (Mt 9 ^s; solid contours), for Mbh = 2.4M© and 
Mms = 40M Q (13.8-Mq CO star). The dashed aslant line shows 
Ei so ,5l = -Bt wi 5i (i?i S o) as is denned by equation (14). Note that 
the possible value of Ei so that leads to Mbh = 2.4Mq is a func- 
tion of M; so (see the main text). Three cases are shown for il- 
lustration, assuming the temporally constant wind injection: If 
r» 1, the required value for Ei so is independent from M so , and 
it is required that Ei so 51 ~ 70 (otherwise the outward shock- 
wave cannot be launched at Mbh = 2.4Mq.) If T w ~ 1, it is 
required that the product -Ej so Mi S o should exceed a certain value, 
i.e., -BisoMso ~ 1 for the outward shock wave being launched at 
M BH = 2.4M Q (e.g., Ei so ,5i ~ 10 if Mso,© = 0.1, and E iso ,5i ~ 1 
if Mso,© ~ 1). 



5.7E ; 

24E; 



1/6 



so, 5 

1/4 



p-1/3 1/12 



-3/4 



iso,52 L sw,8 



for t <t v 
for t > t v 



(13) 



This expression explicitly utilizes the assumption that Ei so 
is constant in time for t ^ t w and zero afterward. 

Since Mbh increases monotonically with time, the den- 
sity structure of the stellar mantle is specified for given 
Mbh (Fig. 2). Therefore, we can estimate the enclosed mass 
within a sphere in which the post-shock materials attain a 
certain temperature (Fig. 2; see also figure 7 of Maeda & 
Nomoto 20 03). 56 Ni is synthesized in the region (Mt 9 ^s) 
where T 9 > 5. Figure 4 shows M Tg ^ 5 [= M( 56 Ni) /fa; the 
isotropic equivalent value for M( 56 Ni)] for Mbh = 2.4M Q 
and Mms = 4OM . 

Mt 9 ^s depends not only on i?iso and E\ so , but also on 
Mbh- On the other hand, as concluded in §3.1, there is a 
condition in terms of £iso and Mi ao , to initiate the explosion 
at given Mbh: This indirectly relates Mi so to Mbh- In order 
to obtain Mbh = 2.4M and M ms = 4OM , it is required 
that -Eiso.siMiso,© > 1 (for T w ~ 1) or that S iBO , 5 i > 70 
(for F w ^> 1). Vertical dotted lines in Figure 4 explicitly 
describe the requirement: (1) for r» 1, -Eiso, 51 = 70, (2) 
for T ~ 1 and Mso,© = 0.1, S is o,5i = 10, and (3) for T ~ 1 
and Mso,© = 1, -Ei so ,5i = 1- 

Interestingly, the dependence of M( 56 Ni) on £^0 and 
-Eiso (Fig. 4) changes the behavior at a specific aslant line in 
the -Eiso — Eiso plane. The line is obtained by equalizing the 
two RHS's of equations (13) and defined as 
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(f) 



j-l/2p 2/3 
P 6 iiiso,51 



(14) 



with Tg = 5, i.e., the temperature necessary for the produc- 
tion of 56 Ni. For given Ei BO , M( 56 Ni) is larger for larger Ei BO 
as long as E iBO ^ E tvl , but saturated for E iBO > E tvr . 

The dividing energy injection rate, E tv ,(Ei BO ), appears, 
because of the different behavior of the post-shock tempera- 
ture before and after t = t w (Maeda 2004). The mass of 56 Ni 
is determined by the radius of the outward shock wave when 
the temperature drops down to Tg — 5. If Ei so ^ E t „ (for 
given £"iao), then the temperature behind the shock wave 
drops down to Tg — 5 during the phase when the energy 
source is still active (i.e., t < t w ). Thus, larger Ei BO corre- 
sponds to the larger amount of energy contained behind the 
shock wave, leading to the larger radius of the shock wave 
when the condition Tg = 5 is satisfied (e.g., see Fig. 2). On 
the other hand, the condition Ei BO > Et„ corresponds to 
the case in which all the energy from the central source is 
liberated in short time scale such that the temperature is 
still high (Tg > 5) at t = t w . The following evolution of the 
shock wave, which determines M Tg ^5, is controlled by i5i so , 
not by Ei so . 

It has been suggested that M( 56 Ni) is larger for larger 
Ei BO (e.g., Nakamura et al. 2001b), but this holds true only if 
the energy is generated and liberated almost promptly (i.e., 
t„ — > and Ei B o — > oo). For example, if Ei so j,i = 1, then 
M( 56 Ni) does not depend on the total energy Ei BO (or E w ) 
(as long as £ iso , 6 i > 0.2) (Fig. 4). 



3.2.2 56 Ni in the wind/ jet 
Mass of the wind is 

which must be smaller than the mass of the accreted materi- 
als, i.e., M w < Mbh — 1.4Mq. Composition of the materials 
within the disk wind is largely uncertain, because it depends 
on the thermal history of the wind in the vicinity of the cen- 
tral remnant (i.e., well below i? w ). The deep understanding 
of the composition of the wind material requires numerical 
calculations or analytic investigation of the innermost part 
of the collapsing star (e.g., MacFadyen & Woosley 1999; 
Pruet, Thompson, & Hoffman 2004; Nagataki et al. 2007). 

However, it is still possible to make a rough estimate, 
in order to see in what regions in the Ei so — Mi so plane 
the disk wind is important. Assuming that the energy of 
the wind is initially thermal energy-dominated near the 
central remnant, the typical entropy of the wind material 
s = S/(fe/TOu) (where &b and m u are the Boltzmann con- 
stant and atomic unit mass, respectively) is written by Ei BO 
and Miso as follows: 



22(i5i SOi 5i Mj so ,0) 



,3/8 



11 - C/ iso,52 



for r v 
for r„ 



~ 1 , 
> 1 . 



(16) 



We find that s > 1 in the parameter range of our in- 
terest; s can be as small as about unity, only when either 
i5i so Mj so or Ei BO is very small, but in such a case the whole 
star collapses onto the central remnant (§3.1). Thus, when 
Ei B o and Mi so are large enough to result in a supernova 
explosion, the wind material should experience the strong 



a-rich freezeout, leaving mainly 4 He, not 56 Ni. The strong 
a-rich freezeout is consistent with numerical calculations of 
the collapsing star (Nagataki et al. 2007). In what follows, 
we take M( 56 Ni) = 0.2M W in the wind which is typical for 
the strong a-rich freezeout (e.g., Pruet et al. 2004; but see 
§5 for further discussion on the uncertainty). 



3.3 Geometry of the ejecta 

The geometry of the bulk supernova materials provide a use- 
ful constraint on the model, since recently more and more 
observational data have been available to address the ejecta 
geometry (see Wang & Wheeler 2008 for a review of spec- 
tropolarimetry; for recent spectroscopy, see Maeda et al. 
2008; Modjaz et al. 2008). Here we give a rough estimate 
of how the geometry depends on the wind parameters. Note 
that the geometry of the bulk supernova materials is differ- 
ent from the geometry of the wind at the injection measured 
by /n; Even if the wind is initially collimated (/n < 1), the 
bulk expansion of the stellar mantle, as is induced by the 
wind, can be less collimated or even quasi-spherical. 

The shock breakout time (t SD ) of the wind/jet is esti- 
mated by 



tsb 



^shock 



-Rco/fshock , 



2-Elso^sb 



(M c 



Mb 



1/2 



(17) 
(18) 



Here Rco and Mco are the radius and the mass of the CO 
star. For M ma = 40M Q , these values are Rco ~ 10 1() cm 
and Mco ~ 13.8M Q . 

As the wind/jet pushes the stellar mantle, the wind/jet 
loses its energy by depositing its energy into the surround- 
ings. If the energy injection is terminated before it breaks 
through the progenitor surface (Rco), a large fraction of the 
energy of the wind/jet are transferred to the stellar mantle. 
This results in a quasi-spherical explosion of the bulk of the 
stellar mantle, even if the wind/jet is initially collimated. 
This happens if the following condition is satisfied. 



Eis 



< tsb 



R 2 co(Mco 



2E is . 



M BH ) \ 



1/3 



(19) 



On the other hand, if £ w > t eo , then a large amount of mate- 
rials are ejected toward the jet direction as compared to the 
equatorial direction, resulting in a strongly jetted explosion. 



4 COMPARISON TO OBSERVATIONS 

The expected SN properties can be expressed by i5i so and 
Mi so and the temporal evolution of these, once the other 
parameters are specified (M ms , -Eiso or E w , and fa). For 
comparison to the observation presented below, we examine 
the simplest case, in which _Bi S o and M so are constant in 
time for t < t w and zero afterward (i.e., Ei so = Ei so tw). This 
is a situation examined in most of the previous numerical 
studies, except for MacFadyen et al. (2001) and Maeda & 
Nomoto (2003) who examined the case where the energy 
injected by the jet (wind) is connected to the accretion rate 
to the central remnant. 

First, the momentum balance determines when the out- 
ward shock wave is launched, yielding Mbh at this time as 
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a function of _Z5i so and Mi so (§3.1; Figure 3). Then, the anal- 
ysis presented in §3.2.1 gives us the mass of 56 Ni, as this is 
given as a function of Mbh, -Eiso, and Ei so (Figure 4). The 
mass of 56 Ni in the disk can be roughly evaluated using the 
result of §3.2.2 (equation 15). Finally, the typical geometri- 
cal feature can be derived as a function of Ei SO , Ei so , Mbh 
(§3.3; equation 19). 

Using the expressions derived in §3, we characterize 
properties of a SN explosion as a function of Ei so and Mi so 
for a given progenitor model (M ms ). Other parameters are 
fn, E w (or t w ), but these can be set without large ambiguity 
by observations for some GRB-SNe of special interest. 



4.1 SN 1998bw: the origin of j6 Ni 

SN Ic 1998bw is a prototypical SN associated with a GRB. 
For SN 1998bw, intensive observational data in the optical 
wavelength are available: modeling these observations (by 
one-dimensional radiation transfer calculations; Iwamoto et 
al. 1998; Woosley et al. 1999; Nakamura et al. 2001a) sug- 
gests that .Ei so ~ 30 and M ms ~ 40M©. Adding to this, there 
exists intensive study for this object using multi-dimensional 
radiation transfer calculations (Maeda 2006a; Maeda et al. 
2006bc; Tanaka et al. 2007). The study suggests that the 
intrinsic explosion energy is smaller than the isotropic value 
by a factor of about 3, inferring that fn ~ 0.3 is a good 
approximation. 

Figure 5a shows the result for i?i SO ,5i = 30, fn = 0.3, 
and M ms = 40M Q . Here, we try to constrain the properties 
of the wind generated by the activity of the central engine 
of SN 1998bw, i.e., Ei so and Mi so (or equivalently, E w and 
M w ). For SN 1998bw, observations constrain three quanti- 
ties, the ejecta mass, M( 56 Ni), and the shape of the ejecta: 
(1) The ejecta contain a large amount of the CO-core mate- 
rials, so that Mbh < WM@, as inferred by the optical light 
curve and spectra (§3.1). (2) M( 56 Ni) ~ O.4M to explain 
its peak luminosity (§3.2). (3) The ejecta are suggested to be 
aspherical, but still a large amount of materials are ejected 
into the equatorial direction as inferred especially by spectra 
at ~ 1 year since the explosion (Patat et al. 2001; Mazzali 
et al. 2001; Maeda et al. 2002). This indicates that the ex- 
panding supernova ejecta are quasi-spherical, rather than 
extremely bipolar (§3.3) 

The three conditions are simultaneously satisfied if 
Ei S o,5iM iso , > 1 and £ iso , 5 i > 20 (for T w ~ 1), or 
£iso,si > 70 (for T w 3> 1) (region A in Figure 5b). In these 
cases, 56 Ni is predominantly produced at the shocked stel- 
lar mantles. The wind origin for 56 Ni is disfavored for SN 
1998bw. The wind contribution exceeds the shocked man- 
tle contribution only if Ei SO ,5i < 5 and Mi so ,© > 0.2 (i.e., 
region C of Figure 5b). However, this combination of the 
parameters results in t v /t so > 1 and an essentially bipolar 
explosion, as is inconsistent with the observation. 

In short, the wind should be massive and long-lived (i.e., 
large Mj ao and large £ w ; see equation 15), in order to pro- 
duce a large amount of 56 Ni within the wind (as large as 
~ 0.4Mq). However, such an explosion with a long-lived 
energy injection results in the extremely bipolar explosion, 
since the outward shock wave can reach the stellar sur- 
face before the energy injection is terminated: This argues 
against the wind-origin of 56 Ni in SN 1998bw. 



4.2 Diversity of GRB-SNe 

Figure 5b shows the expected characteristics of GRB-SNe for 
Ei BO = 30, fn — 0.3, and M ms = 4OM . A variety of features 
are predicted for the wind-driven supernovae depending on 
E iso and M iso . 

If £ is „,5iM is o,0 > 1 (r w ~ 1) or £ iso , 51 > 70 (r w > 1), 
then the resulting SNe should be similar to SN 1998bw in 
the ejected mass. The features are further divided as follows 
(regions A, B, and C in Figure 5b): 

(A) SN 1998bw-like - GRB-SNe 1998bw, 2003dh, 
and 20031w: In the parameter region A, the wind-driven 
supernova is similar to GRB-SN 1998bw, in virtually all the 
observed characteristics. The luminosity is similar to that of 
SN 1998bw, since M( 56 Ni) ~ 0.3 - O.6M . 

This category can account for GRB-SNe 1998bw, 2003dh, 
and 20031w. The region is relatively large in terms of Ei so 
and Mi so , and the features of supernovae are insensitive to 
-Eiso and Mi so ; This can explain why these three GRB-SNe 
are similar in the optical properties. 

(B) Sub-luminous SN 1998bw-like - GRBs 040924 
and 041006?: In the parameter region B, the expected su- 
pernovae are similar to SN 1998bw in the ejected mass and 
the geometry, but the difference is seen in M( 56 Ni). This 
indicates that the supernovae has similar optical properties 
with SN 1998bw, with the diversity in the luminosity, cover- 
ing ~ 0.3 — 1 times that of SN 1998bw. This diversity arises 
because, in this parameter region, Ei so is smaller than Et w , 
and thus M( 56 Ni) is dependent on E iB o (§3.2.1) [note that 
-Eso is larger than E tv , in region A, and thus M( 56 Ni) is 
independent from Ei so unlike region B]. This parameter re- 
gion B can account for sub-luminous SNe found as bumps 
in optical afterglows of some GRBs (e.g., GRBs 040924 and 
041006). 

(C) Bipolar SN 1998bw-like: In the parameter region 

C, a resulting SN is similar to SN 1998bw in its luminosity 
and ejecta mass, but the ejecta are more highly beamed 
than SN 1998bw. We have never directly observed such pe- 
culiar GRB-SNe. If such a GRB-SN is observed in future, it 
will provide a strong evidence of the diverse property of the 
central engine, in terms of Ei so and Mi so - 

Here, 56 Ni is mainly produced within the wind materials. 
The bipolar explosion results from the central energy source 
being active for a long time period (c.f., equation 19; e.g., 
low Eiso smaller than E w ). As such, (1) the mass of 56 Ni 
produced within the shocked stellar mantle is small, and (2) 
the wind should be massive to initiate the explosion (to pro- 
vide the sufficiently large momentum in the non-relativistic 
regime). As a result, the wind contribution dominates over 
the shocked stellar mantle in the production of 56 Ni for the 
(extremely) bipolar explosion case. 

If £iso,5iMiso, < 1 (r w ~ 1) or E iso , 5 i < 70 (r w > 1), 
then the resulting SNe should intrinsically different from SN 
1998bw in the ejected mass (§3.1; regions D and E in Figure 
5b). At the low end of Ei so or Mi so , a supernova can never 
be triggered. 

(D) Bipolar SN with the ejected mass < IMq — 
GRBs 060505 and 060614?: In the parameter region 

D, only a small amount of materials (< IMq) near the sur- 
face of the CO layer are ejected. The supernova ejecta are 
essentially bipolar in this case. 56 Ni is mainly originated in 
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Eiso, 51 Eiso, 51 



Figure 5. (a) Estimate of M( 56 Ni) in the E iso — M iso plane. The progenitor is the 13.8M Q CO star (M ms = 40M Q ), taken from 
Nomoto & Hashimoto (1988). The other two parameters are set as Ej S o,5l = 30 and fa = 0.3, representing the bright GRB-SN 1998bw. 
M( 56 Ni)/MQ is shown for the shocked stellar mantle (red contours) and the wind (blue contours). Mbh/ Mq is shown by the black-dotted 
contours. The shaded region denoted by "Forbidden" is excluded by the requirement that the ejected mass must not exceed the available 
accreted mass budget. The dark shaded region denoted by "Collapse" corresponds to whole collapse of the CO star. The green curve 
shows the line where t w ~ i s b- (b) Expected characteristics of GRB-SNe for Mms = 40Mq. In each region, "L" denote a rough estimate 
of the SN peak luminosity normalized by that of SN 1998bw. 



the wind. Thus, M( 56 Ni) oc M w /£w, ranging from M( 56 Ni) 
< 1CP 5 M0 to ~ O.3M . The expected luminosity is thus 
diverse, depending on E[ so and M[ so . This may correspond 
to non-detection of supernova features in GRBs 060505 and 
060614 (which may also be explained by region E below). 

(E) No supernova - GRBs 060505 and 060614?: In 
the parameter region E, the wind injection can never set the 
outward shock wave. This case corresponds to whole collapse 
of the progenitor CO star without a SN. 



5 CONCLUSIONS AND DISCUSSION 

In this paper, we discussed theoretically expected charac- 
teristics of a supernova driven by the wind/jet. We found 
that the resulting supernova features can be categorized as 
a function of Ei so and Mi sa - Thus, it is possible to constrain 
the nature of the central engine of GRBs by observations of 
associated supernovae. 

Results of the past numerical studies can be understood 
as a limiting case. For example, Nakamura et al. (2001b) 
showed that M( 56 Ni) is larger for larger Eiso. In this paper, 
we have clarified that it holds true only for the limiting case 
where t w — -> and thus Ei so — * oo. For an another exam- 
ple, Tominaga et al. (2007b) showed that Ei so determines 
M( 56 Ni). In this paper, we have shown that this behavior 
appears as a limiting case where F ^> 1, and the dependence 
is different when r ~ 1. 

For SN 1998bw, we find that observations are repro- 
duced only if S iso ,5iMi SOiQ > 1 and E iso , 51 > 20 (if r w ~ 1), 
or Siso.si > 70 (if T w > 1) (region A in Fig. 5 b). We favor 
the shocked stellar mantle as the main site of the production 
of 56 Ni. 

Furthermore, the observed diversity of supernovae asso- 
ciated (or not associated) with GRBs can be accounted for 
by the diverse properties of the wind from the central rem- 
nant (iJiso and Mi ao ). It is shown that the different wind 



properties can potentially explain the diversity of super- 
novae associated with GRBs. 

Further diversity can arise from different progenitors 
(M ms ). Irrespective of M ms , the expected features of the 
wind-driven supernovae can be categorized as we did for 
Mm S = 40M© (§4), i.e., regions A - E. The positions of the 
boundaries between different regions, as well as maximum 
M( 56 Ni) at the high end of E iso and M iso , are dependent 
on M ms . For M ms = 25Mq, the boundary between regions 
A and D is located at -Eiso.siMiso,© ~ 0.05 (for F w ~ 1) or 
-Eiso.si ~ 20 (for T w > 1), smaller than for M ms = 40M Q by 
about one order of magnitude (Fig. 3). Using fa = 0.3 and 
E iso ,Bi = 30, maximum M ( 56 Ni) is ~ 0.24M Q even if M BH = 
1.4M©. Thus, M ms = 25M never yield bright SN 1998bw- 
like SNe. We therefore confirmed that the prototypical GRB- 
SN 1998bw and the similar SNe 2003dh and 20031w must 
have originated from a massive progenitor (Mms ~ 40M©) 
from the nucleosynthesis argument. 

Also, the diversity arising from M ms may be important 
to explain the sub-luminous supernovae possibly associated 
with GRBs 040924 and 041006. These can be explained 
by M ma ~ 4OM with smaller E iso and M iso than for SN 
1998bw. Alternatively, it is also possible that the central en- 
gine provides i5i so and Mi so similar to SN 1998bw, but the 
progenitor mass is smaller than SN 1998bw (i.e., region A 
for M ms = 25M©). It is thus important to derive M ms for 
these sub-luminous cases, by spectroscopic and light curve 
modeling, to distinguish these possibilities. 

On the other hand, the diversity in terms of Mms is 
not important in interpreting no detection of supernovae 
associated with GRBs 060505 and 060614. If they are the 
outcome of the core collapse of a massive star, the properties 
of the wind/jet should fall in the low end of J5 W and M w 
(region E), or region D with Mi so ,0 < 10~ 3 . If this is not 
the case, we should have detected associated supernovae. 
Note that Tominaga et al. (2007b) showed that the wind/jet- 
driven explosion with small i5i S o "can" account for the non- 
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detection of supernovae in these GRBs, while this study has 
shown that E[ so and/or M[ so "must" be small in the central 
engine of these GRBs. These GRBs highlight our suggestion 
that we can constrain the properties of the central engine 
by observations (non-detection, in this case) of associated 
supernovae. 

The present analysis is based on one-dimensional calcu- 
lations. We predict that the different categories (A-E) are 
divided rather sharply in terms of Mi so and -Ei so , while in 
reality this may well be smoothed by the jet-accretion in- 
teraction (Maeda & Nomoto 2003; Tominaga et al. 2007). 
Although the jet-accretion interaction should create a non- 
radial flow, key conclusions in the present paper are not 
defeated. For example, our model can explain the qualita- 
tive behaviors found in the past numerical study, and can 
reach at least rough quantitative agreement for the impor- 
tant quantity like M( Ni) (for limited parameter space as 
investigated by numerical study; see also Maeda & Nomoto 
2003). 

The largest uncertainty involved in the present analy- 
sis is in the treatment of nucleosynthesis within the wind 
material, specifically, the assumption that 20% of the wind 
materials become 56 Ni. Surman, McLaughlin, & Hix (2006) 
concluded that the 56 Ni mass fraction could be as large as 
~ 50% for the wind/jet with the entropy in the range be- 
tween 10 and 30. This does not affect our conclusion for 
r w 3> 1, as the mass of the wind materials is anyway small 
in this case. The only conclusion possibly affected by this 
uncertainty is the supernova feature in the region B. From 
equation 16, we see that the entropy of the wind materials 
falls into the range s = 10 — 30 if -Eiao,5iMi SOj ~ 1. This line 
crosses the vicinity of region B, and lead to the contribution 
of the wind materials for the 56 Ni production comparable to 
that of the shocked mantle in this region. As such, M( 56 Ni) 
is this region may have a significant contribution from the 
wind materials. 
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